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Abstract 

A three- and five-component nonlinear Sclirodinger-type models, which describe 
spinor Bose-Einstein condensates (EEC's) with hyperfine structures F = 1 and 
F = 2 respectively, are studied. These models for particular values of the cou- 
pling constants are integrable by the inverse scattering method. They are related 
to symmetric spaces of BD.I-type ~ S0(2r + 1)/S0(2) x S0(2r — 1) for r = 2 and 
r = 3. Using conveniently modified Zakharov-Shabat dressing procedure we obtain 
different types of soliton solutions. 
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1 Introduction 

The dynamics of spinor BECs is described by a three-component Gross-Pitaevskii 
(GP) system of equations. In the one-dimensional approximation the GP sys- 
tem goes into the following multicomponent nonlinear Schrodinger (MNLS) 
equation in ID x-space [1]: 

idt<!>, + dl<^i + 2(|$i|2 + 2|<l>on$i + 2^1$^ = 0, 
idt^o + dl^o + 2(|$_i|' + |$o|' + l^iH^o + 2$*$i<l>_i = 0, (1) 
+ 9^<l>_i + 2(|<l>_i|2 + 2|<l>on$-i + 2<l>t$^ = 0. 

We consider BECs of alkali atoms in the F = 1 hyperfine state, elongated in x 
direction and confined in the transverse directions y, z by purely optical means. 
Thus the assembly of atoms in the F = \ hyperfine state can be described 
by a normalized spinor wave vector t) = t), $0(2;, i), t))"^ 

whose components are labeled by the values of mp = 1,0,-1. The above 
model is integrable by means of the inverse scattering transform method [1]. 
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It also allows an exact description of the dynamics and interaction of bright 
solitons with spin degrees of freedom. Matter-wave solitons are expected to 
be useful in atom laser, atom interferometry and coherent atom transport. 
It could contribute to the realization of quantum information processing or 
computation, as a part of new field of atom optics. 

Lax pairs and geometrical interpretation of the MNLS models related to sym- 
metric spaces (including the model (1)) are given in [2]. Darboux transforma- 
tion for this special integrable model is developed in [3]. In [4] the authors 
study soliton solutions for the multicomponent Gross-Pitaevskii equation for 
F — 2 spinor condensate by two different methods assuming single-mode 
amplitudes and by generalizing Hirota's direct method for multicomponent 
systems. They point out the importance of integrable cases, which take place 
for particular choices of the coupling constants. 

The aim of present paper is to show that both systems mentioned above are in- 
tegrable by the inverse scattering method and are related to symmetric spaces 
[5] BD.I-type: ~ S0(2r + 1)/S0(2) x S0(2r- 1) with r = 2 and r = 3 re- 
spectively. In Section 2 we formulate the Lax representations for the models. 
Section 3 is devoted to the F — 2 BEG model. In Section 4 we construct 
the fundamental analytic solutions of the corresponding Lax operator L and 
reduce the inverse scattering problem (ISP) for L to a Riemann-Hilbert prob- 
lem (RHP). Using the special properties of the BD.I symmetric spaces we 
also obtain the minimal sets of scattering data Tj each of which allow one 
to reconstruct both the scattering matrix T(A) and the corresponding poten- 
tial Q{x,t). This allows us to derive in Section 5 their soliton solutions using 
suitable modification of the Zakharov-Shabat dressing method, proposed in 



2 Multicomponent nonlinear Schrodinger equations for BD.I. se- 
ries of symmetric spaces 

MNLS equations for the BD.I. series of symmetric spaces (algebras of the 
type so{2r + 1) and J dual to ei) have the Lax representation [L, M] = as 
follows 



[13,6]. 



Lip{x, t, A) 
MV^(x,i, A) 



id^ij + (Q(,T, t) - \J)^j{x, t, A) = 0. (2) 
idtip + {Vo{x, t) + XVi{x, t) - \^J)i){x, t, A) = 0, (3) 

Q{x,t), Vo{x,t)^iadf^ + ^[adfQ,Q{x,t)]. (4) 
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where 



f ^ 



Q = p Soq 

0]fsQ 

The 2r — 1- vectors q and p have the form 



J = diag(l,0, ...0, -1). 



(5) 



g = (?2, ■■■,qr, Qr+l, qr+2, q2rT , P = (P2, ■ ■ ■ ,Pr,Pr+uPr+2, ■ ■ ■ , P2rV , 

while the matrix sq represents the metric involved in the definition of so(2r— 1), 
therefore it is related to the metric So associated with so(2r+l) in the following 
manner 



2r+l 

'S'o = XI i~^)'^^^^k,2r+2-k 
k=l 



^0 1^ 



-So 

1 



(Ekn) 



(6) 



Next we will use 

El — (-B±(ei-e2)) ■ ■ ■ ) E±^ei-er)j -E'±ei , -E'±(ei+e,.), ■ ■ ■ , -E'±(ei+e2)) ) 

We will use also the "scalar product" 

{q- E^) = ^{qk{x, t)Ee.,-ek + q2r-k+2{x, t)Ee-,+ek) + qr+l{x, t)E, 



(7) 



ei ■ 



k=2 



Then the generic form of the potentials Q{x, t) related to these type of sym- 
metric spaces is 



Q{x, t) = {q{x, t) ■ Et) + t) ■ Ef), 



(8) 



where E^ are the Weyl generators of the corresponding Lie algebra (see [5] 
for details) and Ajj*" is the set of all positive roots of so(2r + 1) such that 
(a, ei) = 1. In fact A^" = {ei, ei ± e^, k = 2, . . . , r}. 

In terms of these notations the generic MNLS type equations connected to 
BD.I. acquire the form 



iqt + qxx + '2{q,p)q- {q, soq)sop = o, 

iPt - Pxx - '^{q,P)p- (P, Sop)soq^ 0, 



(9) 



In the case of r = 2 if we impose the reduction pk = ql and introduce the new 
variables $i = q2, $o = qs/V^,, $-i = qi then we reproduce the equations 
(!)• 
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3 F=2 spinor Bose-Einstein condensate, integrable case 



Let us introduce Hamiltonian for MNLS equations (9) with p — e^, e — ±1 

i^MNLS = y dx (^{dj^q, dj^q*) - e{q, (ff + e(g, s^q) {q*, Soq*)^ , (10) 
Define the number density and the singlet-pair amphtude by [7,8,4] 

n =(#,#*)= Yl *a*a, e = (*,So*). (11) 

a=-2...2 

where $2 = q2, $1 = qs, $0 = Qa, ^-1 — Q5,^-2 — Qe- Then the singlet-pair 
amplitude take the form [7,8,4] 

e = 2$2^-2 - 2$i$-i + $0- (12) 

The physical meaning of 6 is a measure of formation of spin-singlet " pairs" of 
bosons. The assembly of atoms in the F — 2 hyperfine state can be described 
by a normalized spinor wave vector 

*(X, t) = ($2(X, t), $i(x, t), ^oix, t), $_i(x, t), $_2(X, t)f, (13) 

whose components are labeled by the values ofm^ = 2,l,0,— 1, —2. Here the 
energy functional within mean-field theory [9,10,7,8,4] is defined by 



The coupling constants Ci are real and can be expressed in terms of a transverse 
confinement radius and a linear combination of the s-wave scattering lengths 
of atoms [1,11,12] and f describe spin densities [4]. Choosing C2 = 0, C4 = 1 
and Co = —2 we obtain integrable by the inverse scattering method model 
with the Hamiltonian. We set for simplicity h = 1, 2m — 1 without any loss 
of generality. The evolution equation is described by the multi-component 
Gross-Pitaevskii equation in one dimension [4] 



.9$ _ 6Egp[^] 
Then we have 



(15) 



4 



or in explicit form by components we have 



(16) 



idt^±2 + d^^^±2 = -2e(d, **)$±2 + e(2$2^-2 - 2$i$_i + $g)$; 
i9t$±i + 9,,$±i = -2e(*, **)$±i - e(2$2^-2 - 2$i$_i + ^l)^ 
idt^o + d^^% = -2e(^, **)$±o + e(2$2$-2 - 2$i$_i + $^)$*. 



4 Inverse scattering method and reconstruction of potential from 
minimal scattering data 



Herein we remind some basic features of the inverse scattering theory appro- 
priate for the special case of F = 2 spinor EEC equations. 

Solving the direct and the inverse scattering problem (ISP) for L uses the Jost 
solutions which are defined by, see [16] and the references therein 



lim 0(x, i,A)e 



iXJx 



U, lim ip{x, t, A)e 



IXJx 



U 



(17) 



and the scattering matrix T(A, t) = ip~^4>{x, t, A). Due to the special choice of 
J and to the fact that the Jost solutions and the scattering matrix take values 
in the group SO{2r + 1) we can use the following block- matrix structure of 
T{X,t) 

I + t-T - \ 

nil ~^ 

r(A,i)= b+ T22 -Sob- , (18) 

^ ct b+'^so ml ^ 

where 6^ (A, t) are 2r — 1-component vectors, T22(A) is a 2r — 1 x 2r — 1 block 
and mf (A), cf (A) are scalar functions satisfying cf = 1/2(6"^ • sob^)/mf. 

Important tools for reducing the ISP to a Riemann-Hilbert problem (RHP) 
are the fundamental analytic solution (FAS) x^^i^^t^ Their construction is 
based on the generalized Gauss decomposition of T(A, i) 



X^{x, t, A) = <p{x, t, \)Sfit, A) = ij{x, t, \)Tf{t, A)D±(A). 



(19) 



Here S'j, Tf upper- and lower- block-triangular matrices, while -Dj(A) are 
block-diagonal matrices with the same block structure as T(A,t) above. Skip- 
ping the details we give the explicit expressions of the Gauss factors in terms 
of the matrix elements of T(A, t) 
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S^{t, A) = exp (±(7^(A, t) ■ Et)) , Tf{t, A) = exp (^(^(A, t) ■ E^)) , 

(20) 



^ mi 0^ 



^ 0^ 



D- 



mj 
\^ 1/mt ) y 

where r^(A, t) = b'^/mf, p*(A, t) = h^/mi and 



012 





m+ = T22 



+ ) 



mo = T22 + 



So 



So 



If Q{x,t) evolves according to (1) then the scattering matrix and its elements 
satisfy the following linear evolution equations 



',—r- 

dt 



±X^b^{t,X) = 0, 



dmf 



dmf 

; 

dt 



0, 



(21) 



so the block-diagonal matrices -D^(A) can be considered as generating func- 
tionals of the integrals of motion. The fact that all (2r — 1)^ matrix elements 
of (A) for A G C± generate integrals of motion reflect the superintegrability 
of the model and are due to the degeneracy of the dispersion law of (1). We 
remind that -Dj (A) allow analytic extension for A e C± and that their zeroes 
and poles determine the discrete eigenvalues of L. 

The FAS for real A are linearly related 

X+(x, i, A) = x-{x, t, X)Gj(X, t), GoA\ t) = SJ(X, t)Sj{X, t). (22) 

One can rewrite eq. (22) in an equivalent form for the FAS ^+(x,t, A) — 
X+(x, t, X)e^^'^^ which satisfy also the relation 



lim ^^ix,t,X) = U. 

A— >oo 



(23) 



Then these FAS satisfy 

e{x,t,x) = r{x,t,x)Gj{x,x,t), 



Gj{x,X,t) = e-^''^G^AX,t)e^ 



(24) 

Obviously the sewing function Gj{x, X, t) is uniquely determined by the Gauss 

factors Sf(X,t). In view of eq. (20) we arrive to the following 
Lemma 1. Let the potential Q{x,t) is such that the Lax operator L has no 
discrete eigenvalues. Then as minimal set of scattering data which determines 
uniquely the scattering matrix T{X,t) and the corresponding potential Q{x,t) 
one can consider either one of the sets %i, i — 1,2 



%,^{p^{X,t),fr{Kt), Ae 



%2 = {r+{X,t),f-{X,t), AeM}. 

(25) 
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Proof, i) Prom the fact that T{X, t) e SO{2r + 1) one can derive that 

= 1 + + hp+,soP+){p-,Sop-) (26) 

for A G M. Using the analyticity properties of mf we can recover them from 
eq. (26) using Cauchy-Plemelji formulae. Given Tj and mf one easily recovers 
?r=(A) and c^(A). In order to recover m2 one again uses their analyticity 
properties, only now the problem reduces to a RHP for functions on S'0(2r+1). 
The details will be presented elsewhere. 

ii) Obviously, given Xj one uniquely recovers the sewing function Gj(a;,t, A). 
In order to recover the corresponding potential Q{x, t) one can use the fact 
that the RHP (24) with canonical normalization has unique solution. Given 
that solution x^(x, A) one recovers Q{x, t) via the formula 

Q{x, t) = hm A (j - x^Jx^ix, t, A)) . (27) 

A— >0O ^ ' 

which is well known. □ 

We impose also the standard reduction, namely assume that(5(x,t) = Q^{x,t), 
or in components pk — q%. As a consequence we have p'{\ t) — p^'*{X, t) and 

T-{X,t)^T+'*{X,t). 



5 Dressing method and soliton solutions 

The main goal of the dressing method [17,18,19,20,21] is, starting from a 
known solutions Xoi^^ '^^ ^) of -^^o(^) with potential Q{o){x. t) to construct new 
singular solutions xti^i t, A) of L with a potential Q(i){x, t) with two additinal 
singularities located at prescribed positions Af ; the reduction p = ensures 
that Xi = {Xi )*. It is related to the regular one by a dressing factor u{x, t, A) 

xf{x,t,X)^u{x,X)x^{x,t,X)uZ\X). u4X)^ lim u{x,X) (28) 

a;— >— oo 

Note that m_(A) is a block-diagonal matrix. The dressing factor u{x,X) must 
satisfy the equation 

idxU + (5(1) (x)u - uQ{Q) {x) - A[ J, u{x, A)] = 0, (29) 

and the normalization condition limA^oo u{x, X) = U. Besides xf (x, A), i = 0, 1 
and u{x,X) must belong to the corresponding Lie group S0{2r + 1,C); in 
addition ^(a;. A) by construction has poles and zeroes at Xf. 
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The construction of u{x, A) is based on an appropriate anzats specifying ex- 
plicitly the form of its A-dependence [6,20] and the references therein. 

u{x,X) = iL + {c{X)-l)P{x,t) + ^-^ - P{x,t), P = S^'P'^So, (30) 

where P{x, t) and P{x, t) are projectors whose rank s can not exceed r and 
which satisfy PP{x, t) — 0. Given a set of s linearly independent polarization 
vectors jn^) spanning the corresponding eigensubspase of L one can define 



P{x,t) = J2 \na{^:t))M^\nl{x,t)l Mab{x,t) = (n^(x, i)|n„(x, i)), 

a,b=l 

X — X+ 

\na{x,t)) = Xo{x,t,X'^)\no,a), c(A) = ^ _ {no,a\So\no^b) = 0. 

Taking the limit A ^ 00 in eq. (29) we get that 

Q^i){x,t)-Q^o){x,t) = {X^ -Xt)[J,P{x,t)-P{x,t)]. 

Below we hst the explicit expressions only for the one-sohton solutions. To 
this end we assume Q(o) — and put Xf — n ^iu. As a, result we get 

q'i'\x,t) = -2iu {P,k{x,t) + {-l)'P-,^,,^,{x,t)) , (32) 

where k — 2r + 2 — k. 

Repeating the above procedure N times we can obtain N soliton solutions. 



5.1 The case of rank one solitons 



In this case s = 1 so that the generic (arbitrary r) one-soliton solution reads 



qk = 

1/^ — II 



u — —2fi, z{x,t) = u^x — ut — ^0), (33) 



_ IIj. ^ '-^fc — I 7 ^~*0 — I I ' 

1^0,1 1 A/|noillno2r+il 2|no,ino,2r+i| 



111 II v^2r I 1 2 

^ J- 1 F0,2r+l| \no,k\ *2 Z^fc=2F0,fc| 

SO — T^in— — , Q^fc — —j= =, ^0 — 



'T'O,!! |?T'0,2r+l| 

and Sq — argno,i//i = — argno,2r-+i/A*) 0fc = a'rgno,^. The polarization vectors 
satisfy the following relation 

E 2(-l)'=+^no,,no,^ + (-l)^<,+i - 0. (34) 
fc=i 
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Thus for r = 2 we identify = q2, = qs/V^ and $3 = and we obtain 
the following solutions for the equation (1) 



3»±i = — T-^ (cos 0±i cosh z±i — i sin ©±1 sinh z±ij (35) 

cosh 2z + Aq 

- - 02 - 04 , 02 + 04 1 , "4 

d±i = do + — , 0±i = ^ 2;±i = 2; T 7^ In — , 

\/2iz/a3e-^^(^-''*-^o) , , . , . . , 

$0 = ; 7^5 (cos 03 smh z — i sm 03 cosh z) . (36) 

cosh 2z + Aq 



For r = 3 we identify $2 = ?2, '^'i = qs, ^0 = ?4, $-1 = q5 and $_2 = Q'e, so 
that the one-soliton solution for equation (16) reads 



2iz/V«^e-^'^(^-''*-''±2) . , , . , . , , 

$±2 = cosh 2;g + ^^^^ ^"^^ ~ * ^"^^^ *^ ^ 

2it/V5ia^e-^/-("-^*-^^0 , . 

$±1 = — — —2 (cos 0±i smh z±i - I sm 0±i cosh z±i) (38) 

cosh 22; + Aq 

02 -06 , 02 + 06 _ 1 1 «6 

C>±2 = do + , 0±2 = ^ ^±2 = 2;+-ln— , 

ZfJ, Z I OL2 

X r _ 03 - 05 , 03 + 05 ^ 1 1 "5 

2/1 2 2 a-i 

2iz/a4e~''^(^~''*~'^°) 

$0 — 79 — (cos 04 cosh z — % sin 04 sinh z\ . (39) 

cosh 2z + Aq 

Choosing appropriately the polarization vectors |n) we are able to reproduce 
the soliton solutions obtained by Wadati et al. both for F = 1 and F — 2 
EEC. 



5.2 The case of rank two solitons 



Here s = 2 and we have two linearly independent polarization vectors \na), 
a — 1,2. Prom eq. (31) we get 

P{x,t) = [\n^{x,t))M22{n{{x,t)\ - \n2{x,t))M^2{n\{x,t)\ 

- \n^{x,t))M2i{nl{x,t)\ + \n2{x,t))M^^{nl{x,t)\) , (^0) 
detM(x,i) =MnM22 -M12M21, Mab{x,t) = {nl{x,t)\nhix,t)) , 

The corresponding expressions for the rank 2 soliton solution are obtained by 
inserting eq. (40) into (32) and are rather involved. We remark here that the 
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reduction = Q may not be sufficient to ensure that det M is positive for 
all X and t. so for certain choices of \na) we may have singular solitons. These 
and other properties of the rank 2 soliton solutions will be analyzed elsewhere. 



6 Conclusions and discussion 

The main result of the present paper is that a special version of the model 
describing F = 2 spinor Bose-Einstein condensate is integrable by the ISM. 
The corresponding Lax representation is naturally related to the symmetric 
space BD.I. ~ SO(7)/SO(2) x S0(5), see [5]. For a generic hyperfine spin F, 
the dynamics within the mean field theory is described by the 2F + 1 com- 
ponent Gross-Pitaevskii equation in one dimension. If all the spin dependent 
interactions vanish and only intensity interaction exists, the multi-component 
Gross-Pitaevskii equation in one dimension is equivalent to the vector nonlin- 
ear Schrodinger equation with 2F -|- 1 components [22] . 

Then equations (9) with the reduction p = eq*,e = ±1 are natural general- 
ization of the vector nonlinear Schrodinger equation, which adequately model 
the spinor Bose-Einstein condensates for values of F = r equal to 1 and 2. We 
expect that for generic F these equations may be useful in describing BECs 
with higher hyperfine structure. 

Here we derived only generic one-sohton solutions. Following the ideas of [23] 
one can classify different types of one-soliton solutions related to different 
possible choices of the rank of P{x,t) and its polarization vectors. One can 
also derive the iV-soliton solutions by either repeating N times the dressing 
with u (see eq. (30), or considering more general dressing factors u with 2N 
zeroes and poles in A. These and other problems will be addressed elsewhere. 
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